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Abstract 

. One of the qualitatively distinct and robust implication of Loop Quantum Gravity (LQG) is the 



underlying discrete structure. In the cosmological context elucidated by Loop Quantum Cosmol- 
ogy (LQC), this is manifested by the Hamiltonian constraint equation being a (partial) difference 



If^ • equation. One obtains an effective Hamiltonian framework by making the continuum approxi- 
mation followed by a WKB approximation. In the large volume regime, these lead to the usual 

m 

I classical Einstein equation which is independent of both the Barbero-Immirzi parameter 7 as well 

■ as h. In this work we present an alternative derivation of the effective Hamiltonian by-passing the 

. continuum approximation step. As a result, the effective Hamiltonian is obtained as a close form 

, expression in 7. These corrections to the Einstein equation can be thought of as corrections due to 



the underlying discrete (spatial) geometry with 7 controlling the size of these corrections. These 



O 

^H^. corrections imply a bound on the rate of change of the volume of the isotropic universe. In most 

cases these are perturbative in nature but for cosmological constant dominated isotropic universe, 

• 1—1 . 

^ ' there are significant deviations. 
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I. INTRODUCTION 



It is widely believed that one needs to construct a quantum theory of gravity. In two of 
the leading approaches, string theory [l| and LQG 1^, the marriage of quantum theory and 
gravity needs some new essential input (s). In the ambitious string theory, one such input 
is extra dimensions still retaining framework of continuum space-time geometry. In LQG, 
the new feature that emerges is the discrete nature of the spatial geometry (as distinct 
from manifold structure). It is also a fact that for most purposes, classical continuum 
geometry and the Einstein equations suffice to an excellent approximation. It is thus a 
natural question: How much deviations from classical general relativity (GR) are predicted 
from the more refined quantum versions? 

In order to address such a question in some quantitative manner, it is necessary to be 
able to identify physical contexts in which the more refined theory can be computed in 
sufficiently explicit terms. A natural context is that of homogeneous cosmologies or even 
simpler, that of homogeneous and isotropic cosmologies. In this context the geometry is 
described in terms of a single variable of a single coordinate - the scale factor - while 
matter can be effectively restricted to scalar field(s) with possibly non-minimal couplings 
(eg dilaton). Even within the cosmological context, there are regimes of small volume (closer 
to the classical singularity) and regimes of large volume where classical GR is expected to be 
an excellent approximation. In this work, we obtain deviations from GR in the large volume 
regime, as implied by Loop Quantum Cosmology 0, 14 1 , as a function of the Barbero-Immirzi 
parameter 7 (and another order 1 dimensionless parameter /ig). 

In the previous works on obtaining an effective Hamiltonian 0, 0, Q , a two step approach 
was taken, namely, a continuum approximation wherein the fundamental difference equation 
is replaced by a second order differential equation (Wheeler-DeWitt equation) followed by 
a WKB approximation of the differential equation. We will refer to the first step as the 
pre-classical approximation [6]. The second order differential equation is independent of the 
Barbero-Immirzi parameter and the Planck constant. One could go beyond the pre-classical 
approximation and obtain higher order (in 7) corrections to the effective Hamiltonian follow- 
ing the same two step approach. In this work, however, we derive an effective Hamiltonian 
as an explicit closed form expression involving 7, by-passing the pre-classical approximation. 

In section |H1 we recall the relevant aspects of isotropic LQG namely, the fundamental 



2 



difference equation, the form of generic solutions, the pre-classical approximation and the 
WKB route to obtaining an effective Hamiltonian constraint. In section IIIIl we make the 
WKB approximation directly at the level of the fundamental equation to arrive at the 
effective Hamiltonian as well as the conditions for validity of the WKB approximation. 
The effective Hamiltonian is obtained as a close form expressions in e := 2/io7, where 
fiQ is a dimensionless parameter of order 1. In section IIVI we discuss some illustrative 
solutions of the effective dynamics for minimally coupled single scalar field. We also discuss, 
how the solutions of the e— dependent equations differ if phenomenological matter (dust, 
radiation and cosmological constant) are incorporated somewhat naively. In the limit e ^ 
these equations and solutions go over to the usual Einsteinian classical solutions. These 
solutions exhibit the situations wherein, the discreteness corrections (codified by e) appear 
as perturbative corrections as well as when the discreteness corrections have qualitatively 
different implications. Finally in|3 we clarify the nature of approximations and domains of 
validity and summarize our results. 



II. SUMMARY OF LQC FRAMEWORK 

n 

We work within the LQC framework as detailed in For our purposes, it may be 
summarized as follows. Classically, we have a phase space with coordinates c,p and the 
Poisson bracket given by {c,p} = where k = IGttG and 7 is the Barbero-Immirzi 
parameter whose precise value does not concern us here. The coordinate c denotes the only 
invariant part of the gravitational gauge connection while p denotes the only invariant part 
of the densitized triad. The loop quantization leads to a kinematical (i.e. before imposition 
of the only surviving Hamiltonian constraint) Hilbert space conveniently described in terms 
of the eigenstates of the densitized triad operator (£| := kH), 

Plfj') = l^^lM . {f^\f^') = ^f^i^' (1) 

Note that the eigenstates of the densitized triad operator are properly normalized even 
though the eigen values span the real line. This is a reflection of the fact that in loop quan- 
tization, only holonomies of connections are well defined operators and not the connections 
themselves. 

As in LQG, the classical Hamiltonian constraint is promoted to a quantum Hamiltonian 
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operator using the volume operator whose action on the triad basis states is given by 

v\^^) 



1 ,2 



7/> 



g,-p.^ l/i) := V,\^^). (2) 

In the isotropic context we have two classes to consider, namely spatially flat and closed 
models. The quantization of the corresponding Hamiltonian operators is given in ^. By 
introducing a parameter r] we can deal both classes together. The values = and rj = 1 
will give the fiat and the closed models respectively. 

In loop quantization of the classical gravitational Hamiltonian (in connection formula- 
tion), one has to express the connections (and their curvatures) in terms of appropriate 
holonomies of the connection. This introduces a dimensionless, real parameter /iq through 
a fiducial length of the loops used in defining the holonomies. Thisparameter is to be de- 
termined relating this symmetry reduced model to the full theory ^. For the purposes of 
this paper, it is simply treated as a parameter of order 1. 

The action of the gravitational Hamiltonian on the triad basis states is then given by 
3 

(e-wr,|^ + 4^^) - {2 + 4fil'y'r])\fi) + e^^°^|/x - 4/io) ) . (3) 

Notice that the Hamiltonian connects states differing in their labels by ±4/io. This is a 
direct consequence of the necessity of using holonomy operators in the quantization of the 
Hamiltonian operator and is responsible for leading to a difference equation below. 

General kinematical states in the triad basis have the form 



where n run over some countable subset of M. 

The Hamiltonian constraint of the classical theory is promoted as a condition to define 
physical states, i.e., 

(^S+^matL)k) = 0. (5) 

In terms of := e~i^Spt and The Hamiltonian constraint © translates into a difference 
equation. 



1 



+ 8/s:7V(^g7^pJ ^m(/i)5^ , V/i eM (6) 

Af, := |/X + /io|^ - |/i-/io|^ , ^matterl/^) := ^m(/^)|/^) • 



Hmifi) is a symbolic eigen value and we have assumed that the matter couples to the 
gravity via the metric component and not through the curvature component. In particular, 
= 0) = 0. A couple of points are note worthy. 

Although yU takes all possible real values, the equation connects the 5^ coefficients only 
in steps of 4/io making it a difference equation for the coefficients. By putting /i := z/ + 
(4/io)n, 72 G Z, z/ G (0,4/io), one can see that one has a continuous infinity of independent 
solutions of the difference equation, labeled by z/, Sn{v) '■= Si.+4^p„. For each u an infinity 
of coefficients, Sn{i')-,^n G Z, are determined by 2 'initial conditions' since the order of the 
difference equation in terms of these coefficients is 2. Coefficients belonging to different v 
are mutually decoupled. Denoting any two (fixed) independent solutions of the difference 
equation by Pri(z/) and cr„(z/), a general solution of the second order difference equation, for 
each V G (0,4yUo) can be expressed as Sniy) = >S'o(z/)p„(z/) + S'i(z/)cr„(z/) where 5*0(2/), 5*1 (z/) 
are arbitrary complex numbers. Linearity of the equation implies that only their ratio is 
relevant and for u = 0, this ratio is fixed 2|- A general solution of the fundamental equation 
(jni) is given by a sum of the 5'„(z/) solutions for a countable set of values of G [0,4yUo). 
Note that there are infinitely many ways of selecting the countable subsets of values of u. 
Since the coefficients and the symbolic eigenvalues Hm{fJ'), both vanish for fi = 0, the 
coefficient sq decouples from all other coefficients (sq of course occurs in the u = sector 
j^). As the are defined for all /i G M, there is no break down of the dynamics at zero 
volume and in this sense, the LQC dynamics is singularity free jlfl| . 

For large values of /i ^ 4/io {n ^ 1), which correspond to large volume, the coefficients 
An become almost constant (up to a common factor of ^/n) and the matter contribution 
is also expected similarly to be almost constant. One then expects the coefficients to vary 
slowly as n is varied which in turn implies that Sn{i^) also vary slowly with n. This suggests 
interpolating these slowly varying sequences of coefficients by slowly varying, sufficiently 
differentiable, /uncizons of the continuous variable p(n) := |7^pn 0. Using Taylor expansion 
of the interpolating function, the difference equation for Sn{i^) then implies a differential 
equation for the interpolating function. This is referred to as a continuum approximation. 
The terms up to (and including) second order derivatives, turn out to be independent of 7 
and the differential equation, truncated to keep only these terms, matches with the usual 
Wheeler-DeWitt equation of quantum cosmology. This will be referred to as a pre-classical 
approximation ^. 
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In the second step, one makes a WKB ansatz for solutions of the differential equation. To 
the leading order in h, one obtains a Hamilton- Jacobi equation for the phase from which the 
effective Hamiltonian is read-off. In 7], such an effective classical Hamiltonian was obtained 
from the pre-classical approximation and it was shown that the classical dynamics implied 
yy the effective Hamiltonian is also singularity free due to a generic occurrence of a bounce 



The derivation of the effective Hamiltonian was based on a WKB approximation for 



the second order differential equation. It is natural to seek to improve the pre-classical 
approximation by including higher derivative terms in the differential equation. 

In this work, however, we will follow a different route. We will make the WKB ansatz at 
the level of the difference equation itself. The slowly varying property of the interpolating 
function will be applied to the amplitude and the phase of the interpolating function. It 
is then very easy to obtain an effective Hamiltonian with nontrivial dependence on the 
Barbero-Immirzi parameter 7 and also the domain of validity of WKB approximation. In 
the limit 7 — > 0, one will recover the earlier results of Q]. 



III. DERIVATION OF EFFECTIVE HAMILTONIAN 

In anticipation of making contact with a classical description, introduce the dimensionful 
variable p{jj,) := ^'jip^ as the continuous variable and an interpolating function ipi^p) via 
ip{p{fi)) := 5^. Correspondingly, define po := |7^p;Uo which provides a convenient scale to 

demarcate different regimes in p. 

3 

Defining A{p) := {^'yip)'^ A^^, q := 4po and replacing by the interpolating function 
ip{p), the fundamental equation © becomes, 

= A{p + q)^{p + q)-(^2 + 9^7]^ A{p)^{p) + A{p - q)i;{p - q) 

+ l^^KHMm (7) 

At this stage one could use 'ipip + q) = YlkLo AT'0fc(p), with ipk{p) denoting the k*^ deriva- 
tive of the function with respect to p, derive an infinite order differential equation, consider 
WKB approximation to extract the o{h) terms which gives a Hamilton- Jacobi equation 
and read-off the effective Hamiltonian. Indeed, this is how we first arrived at the effective 
Hamiltonian, but a much shorter route is available. 
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The equation (|7j) can be thought of as a functional equation determining ipip)- ^ind of 
course it has infinitely many solutions, corresponding to {Sn{i^)},i^ € [0,4/io). So far there 

has been no approximation, we have only used p{fi) instead of fi and ip{p) instead of s^. 

Mp) 

One can always write the complex function ip{p) := C{p)e^ ^ . In the spirit of a contin- 
uum description, we assume that the amplitude and the phase are Taylor expandable and 
write, 



C{p±q) 
5C+ 

5$, 



C{p){5C+±5C. 

n=0 

00 

E 



[2n)\ Cip) 



^ (2n)! 

n=l ^ ' 



$(p±g) := + 5$+ ± 



00 



n=0 



(2n + l)! C(p) 



(8) 
(9) 

(10) 



We have just separated the even and odd number of derivative terms for later convenience. 
Substitution of ^{p ± q) = C{p ± g)e^*(P^'') in © leads to. 







Hr, 



1 6 



(2 + e^r])A{p) + 



1 6 



(B^p, q)6C4p) + B_ip, q)5C.{p)) [cos 
iB.ip,q)SC+ip) + B+ip,q)5C.ip)) (si 



T' 

sm— — 
V ^ 
(5$. 



cos- 



sm- 



+i{B_{p,q)5C+{p) + B+{p,q)5C^{p)) cos— ^ sin 



+ (sin 



h 



cos- 



B±ip,q) := A(p + g) ± A(p - g) 



Thus, one has two equations for the /onr combinations, SC±, 

Now we assume that the amplitude and the phase are slowly varying functions of p over a 
range q i.e. when compared over a range ±g, successive terms of a Taylor expansion about p 
are smaller than the preceeding terms. For example, \^\ -C 1, \ ^^\ ^ 1 etc. and similarly 
for the phase. 



For such slowly varying solutions of ((Zj), we can approximate the equation by keeping 
only the first non-trivial derivatives i.e. 6C+ = 1, 6C- = ^ and = 5$+ = ^-y-. 
To arrive at a Hamilton- Jacobi equation, we identify ^'{p) := which implies that ^ = 
2/io7^P* _ where we have used the definitions e = 2/io7 and £p = nh. This combination 
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of the first derivative of the phase thus has no exphcit dependence on h. 



The terms in the real and imaginary equations can be organized according to powers of 
h. For sine and cosine of 5$+, we have to use the power series expansions and keep only the 
leading powers of h. For sine and cosine of no expansion is needed since there is no h 
dependence. The real and imaginary equations can then be written as, 



= Hm--4^{2 + e'v)A{p)+^ ^ 



(12) 



_ f q'^"{p) \ , B_{p,q) qC'\ , f B4p,q) , qC'\ tan jeK) 

- [ 2heK ) y + B^{p,q) C ] + W{p,q) ^ C ) eK ^'^^ 

On physical grounds, one expects a classical approximation {hP terms), to be valid only 
for scales larger than the quantum geometry scale set by q and therefore we limit to the 
regime p > q. In this regime, the coefficients A{p), B±{p, q) behave as 3|, 

3 3 
A{p) ^ -gv^-o(p"2) 

3 3 

B+ip,q) ~ -^qVp-o{p'2) 
3 1 

B-{p,q) ~ -g^ + o(p"5) (14) 

Noting that q = ^eip and keeping only the leading powers of ip (or h), the real and 
imaginary equations become (WKB approximation), 

= - ^i(2 + ^'v)Aip) + \^ g)cos(eir)] (15) 

_ ( q^^'\p) \ ( B_{p,q) qC'\ tan (eiT) 

V J^\B^{p,q)^ C J eK ^ 

The real equation, ()15p is of o{hP) and is a Hamilton- Jacobi equation for the phase. The 
right hand side, viewed as a function of p, K, it is the effective Hamiltonian constraint. The 
imaginary equation, (|T6|) is of o{h) and is a differential equation for the amplitude, given the 
phase determined by the real equation. The equations for the phase and the amplitude are 
decoupled. For self consistency of the WKB approximation, the solutions have to be slowly 
varying. In general, the solutions will be slowly varying only over some intervals along the 
p-axis and such interval (s) will be the domain of validity of the WKB approximation. One 
can infer the domain of validity as follows. 

Consider the eqn. (fT6|) . By definition of slowly varying phase, the absolute value of the 
first term must be much smaller than 1. Since smallest value of ^'^^j^^ is 1, (the absolute 



value of) the bracket in the second term must be smaller than 1. Since ^ is also small by 



definition of slowly varying amplitude, we must have | < 1. This immediately requires 
p > q = and is obviously true in the regime under consideration. Evidently, for smaller 
volumes, p ^ q, where ^ is larger, we must have > 1 i.e. eK ^ 0. By contrast, for 

larger volumes, p ^ q, larger values of eK{< |) are permitted. Thus, the equation (fTHj) 
serves to identify a region of the classical phase space where the effective classical description 
is valid. The domain of validity of effective Hamiltonian and its relation to the usual general 
relativity (GR) Hamiltonian is discussed further in the remarks below. 
It is convenient to write the real equation in the form. 



= 



+ - 



\^B4p,q) (^^sin^ j +2A{p)v 

{B+{p,q)-2A{p)} 
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+ H„ 



(17) 



The second square bracket in eq. ()17jl is precisely the quantum geometry potential, W^g/n, 
It is unaffected by the implicit inclusion of all terms of the Taylor expansion of ip{p). 
Since the leading terms in — 2A vanish, the quantum geometry potential term is inde- 
pendent of e. But it is also higher order in h. From now on, this term is suppressed. It is 
obvious that in the limit e — > 0, one gets back the expression for the effective Hamiltonian 
constraint obtained from the pre-classical approximation in 

Let us focus on the effective Hamiltonian constraint in the large volume regime. As 
mentioned earlier, the matter is assumed to couple only through metric components and 
thus the matter Hamiltonian has dependence on p but not on K. For the moment, we will 
not be explicit about the matter Hamiltonian and simply view it to be a function p and 
matter degrees of freedom symbolically denoted by (f),p^. In the next section we will be 
more explicit about it. The relevant Hamiltonian constraint then becomes (suppressing the 
quantum geometry potential): 

'A 







H„i{p, 



(18) 



To see what the modified Hamiltonian constraint implies for cosmological space-times, 

we have to obtain and solve the Hamilton's equations. To this end, we identify \p\ = jo^ and 

choose the synchronous time as the evolution parameter (Lapse = 1). It is straight forward 

to obtain the Hamilton's equations of motion as: 
dp ^/p . . ad 



dt 



(19) 



dK _ K dHm{(j),p^,p) ^ 1 
dt 3 dp Ay/p 

d4 _ dHm{p,4>,P<t>) 
dt dp^ 

dHm{p,(j),p<p) 



rsm 



K) +r] 



dt 



(20) 
(21) 
(22) 



The equation pSj) is the modification of the Friedmann equation while (fT^ OUj) lead to 
modified Raychoudhuri equation after eliminating K, p in favour of the scale factor and its 
time derivatives. In particular, (fTT?|l gives d = and leads to cos(e-ft') = y/l — e^d^ 

(e — s> limit fixes the sign of the square root). Thus the Hamiltonian constraint can be 



expressed in terms of a. Furthermore, since the constraint is obviously preserved along 
the solutions of the Hamilton's equations, we can obtain K in terms of H^- Now is it 
is straightforward to construct left hand sides of the Friedmann and the Raychoudhuri 
equation and comparing with usual Einstein equations, read-off the effective density and 
pressure. One gets. 



K 



a? + rj 



K 

2 



8a -^Hrr,. - —Hm - V 



(pefr + 3Pefr) 



a K 



[Hm. — a 



2^0^ \3a 
dHr 



da 



2 U""""^ 



(23) 

(24) 



For future reference, we note that the effective density and pressure can be expressed as ^ 



Pes 

(PefT + 3Pefr) 



8a ^Hm 

8a ^^{^Hrn 



3e'a (Ak 



where, 



(25) 



16k \3a 

As e — i> (7—^0), the right hand sides of the above equations go over to the equations 
of and so do the effective density and effective pressure, with the quantum geometry 
potential terms suppressed. 

The advantage now is that the Raychoudhuri equation is automatically satisfied once 
the Hamilton's equations for the matter hold and the Friedmann equation holds. Thus 



^ Our normalizations are inherited from the normahzations used in Bianchi models, by equating the three 
triad variables \pi\ to the triad variable for the isotropic specialization \p\ which equals The Bianchi 
normalized volume thus equals A Hamiltonian is thus ^ times the Hamiltonian density 
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we have the e-corrected, three coupled (for a single matter degree of freedom), first order 
ordinary differential equations which go over to the classical equations when e — > 0. It is 
now a simple task to compare the solutions with/without e corrections, for the same initial 
conditions consistent with small a and large a. Some solutions are discussed in the next 
section. 

Several remarks are in order: 

(1) It is somewhat surprising that the parameters related to quantum geometry, yUo and 
7 (which have to be non-zero), appear in the effective Hamiltonian which is independent oi 
h i.e., in a c/asszca/ description. 

To trace how this happens, recall that the basic phase space variables in isotropic LQC 
(in the connection formulation) are c,p and the conjugate momentum variable K is related 
to the connection variable c by c = 'jK (for spatially flat model for definiteness) . In loop 
quantizing the Hamiltonian constraint, the classical constraint is to be expressed in terms 
of the holonomies of the connection (equation (34) of [4]). In this process, the parameter 
fiQ enters. The classical constraint so obtained, precisely contains the sin^/xoc which is the 
same as sin^^. The parameter e is thus present already at the classical level (JtP). Within 
a strictly classical context, one can view it as a "regulator" and remove it at will, by taking 
e — s> to recover the classical Einstein Hamiltonian. Within a strictly WKB context also one 
can remove it since the corrections to WKB are positive powers of | ~ o(e). However, from 
the perspective of the exact loop quantization, we cannot take e — 4]. While viewing the 
WKB solutions as an approximation to the exact LQC solutions, we retain the link to the 
exact solutions by keeping e 7^ 0. 

(2) The classical GR Hamiltonian is the expression obtained from the effective Hamil- 
tonian in the limit eK (and large volume). This could be achieved by either taking 
e ^ 0,K fixed, which is a natural classical GR perspective since the classical GR has no 
e or by taking K ^ 0,e fixed which is the appropriate quantum perspective in LQC. This 
opens up the possibility that there could be modifications of classical GR especially when 
the conjugate momentum (K) gets somewhat larger. As noted already, for smaller volumes, 
the domain of WKB approximation restricted by equation (fT^ . requires K to be small and 
the effective Hamiltonian reduces to the GR Hamiltonian. For small K, K ^ —a and cor- 
responds to the extrinsic curvature (a geometrical quantity). For larger values of K, the 
relation between K and the extrinsic curvature is given by equation 
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(3) Since for smaller volumes, the effective Hamiltonian goes over to the GR Hamiltonian, 
the effective density and pressure defined in (12211211) go over to those defined in [3]. Conse- 
quently, the genericness of inflation is insensitive to the e parameter. Furthermore, since 
the quantum geometry potential is unaffected, the genericness of bounce also continues to 



no 



d (in this regime, the quantum geometry potential must be retained in equation (fTTj) ) 



(4) The occurrence of the trigonometric function of K immediately implies two bounds. 
The effective Hamiltonian constraint, ()18p. implies that is that the matter Hamiltonian is 
necessarily bounded: rj < ^'^'^Hm < f] + ^ oi equivalently, < ^a'^Hm — rj < ^. 
Equation on the other hand implies a bound on the expansion rate: \a\ < e^^. 

3 

Consider the bound on the matter Hamiltonian. Momentarily, let us write Hm = 
where the 'density' D may include modifications from inverse triad operator We have 
the two regimes of large and small volumes, two possibilities for spatial curvature, t] and 
two bounds on the matter Hamiltonian to consider. For large volumes, we may identify D 
with the usual energy densities for dust and radiation (late time dominant matter) or with a 
positive cosmological constant. For dust and radiation, the upper bound is satisfied for both 
flat and closed models but the lower bound will be violated for closed models [rj = 1) and 
hence a re-collapse must commence. If D is identified with a positive cosmological constant, 
then the lower bound will be satisfied but the upper bound will be violated for large enough 
volumes irrespective of rj, resulting in re-collapse. For smaller volumes, cosmological constant 
will satisfy the upper bound but will violate the lower bound for a close model implying 
a bounce. In this regime, the dust/radiation will potentially violate the upper bound. 
However, for smaller volumes, D has modifications from the inverse volume operator which 
imply that D goes to a constant ensuring that the upper bound is indeed satisfied. In the 
small volume regime, one also has a contribution from the quantum geometry potential 
which implies that irrespective of r], the lower bound will be violated implying a bounce 
Since the effective constraint excludes evolutions which violate either of the bounds, 
a violation in the large volume regime must encounter a re-collapse while a violation in 
the small volume regime indicates a bounce. The solutions discussed in the next section 
illustrate these implications. 

Consider now the bound on the expansion rate, < a < e~^. Let us trace a possible 
scenario beginning with a bounce for a spatially fiat model. Let us set t = at the bounce 



12 



time and a(0) := cibounce ^ ^p- Near this time, we must include the quantum geometry 
potential i.e. Hm — * Hm '■= Hm + in the equation ()23|) and Hm ~ ha^, where h is 
a constant The bounce is determined by Hm = and gives ~ h^^^^ 11]. Close 

to t ~ 0, the scale factor will behave as, a{t) ^ ctbounce + ott^,(3 > 2 and a is a positive 
constant. Since a = a(3t^~^, the bound will be violated in a finite time. The bounce phase 
will go over to an exponentially expanding phase when the quantum geometry potential 
becomes negligible. This will happen well before saturating the bound on d, even if abounce 
is just a few times £p. Suppose the bounce phase ends at time ti and the universe enters a 
phase of exponential expansion, a{t) ~ e'**. The exponent is determined by the parameters 
and the duration ti of the preceeding bounce phase. Now a{t) = a{t)h which increases 
monotonically. This phase must again end in a finite time (at t2 say) to satisfy the bound. 
Since = ^ and (2(^2) < e-\ the maximum expansion is effectively controlled by a{ti). 
Once again, the exponential inflation can end before saturating the bound and these details 
will depend on the precise matter dynamics (eg the form of scalar field potential). One can 
repeat the scenario for power law accelerated expansion and have a similar finite duration 
phase. If the bound gets saturated at any stage, the acceleration must go over to deceleration 
possibly reaching a re-collapse, contraction and another bounce and so on. If however the 
infiationary phase ends prior to saturating the bound and the universe enters a decelerated 
phase such as radiation domination followed by matter domination, the bound will continue 
to hold. Note that such a scenario is also conceivable without a bound on d, however with a 
bounded d the scenario becomes more plausible. The bound is also specified in terms of the 
discreteness parameter, e, which is expected to be constrained by the full (inhomogeneous) 
loop quantum gravity. 

(5) Finally, let us note that d is the extrinsic curvature of the symmetry adapted hy- 
persurface (and is also the rate of change of the physical volume of the universe) and 
is gauge invariant. The modifications to GR due to the non-zero parameter e are man- 
ifested in the modified coupling (quadratic in matter Hamiltonian) between matter and 
geometry. The GR domain, contained within the domain of validity of WKB, can also be 
given as: \p\ = ^ ^ q = -y' translates into a ^ \f^^P while \eK\ <^ | translates into 
< ^Hjn — Tj ^ The effective Hamiltonian and subsequent analysis being o{h^), is 
insensitive to factor ordering in the Hamiltonian constraint. 
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IV. MODIFIED SOLUTIONS 



Let us consider a few solutions of the modified equations. These are not meant to be 
phenomenologically reahstic solutions, but are to be viewed as indicating modification to GR 
solutions. In particular we consider the cases of a minimally coupled homogeneous scalar 
field, phenomenological matter with a constant equation of state and positive cosmological 
constant viewed as a special case of phenomenological matter. 

A. Minimally coupled massive scalar field 

For simplicity, let us take the matter sector consisting of a scalar field, 0, minimally 
coupled to gravity. Then its usual classical Hamiltonian is given by Hm{(f), p^, p) = 

_3 

p^i'^V{(j)). Its quantization involves two parts - (a) quantization of p 2 operator and (b) 
Bohr quantization of the scalar field itself. For large volume, the inverse volume operator 
in the triad representation just goes over to the classical expression. Its correction terms 
involve higher powers of h (the terms suppressed by inverse powers of p in the coefficients 
B+{Piq)i A{p,q), likewise involve higher powers of K). Thus, to o{h^), these modifications 
are irrelevant in the large volume regime. The Bohr quantization appears not to introduce 
7— corrections to the usual classical matter Hamiltonian. In this paper we will pretend that 
the matter is quantized in the usual Schrodinger quantization. Loop quantization of matter 
with or without non-minimal coupling (to mimic a dilaton) will be presented elsewhere. 

The solutions are obtained numerically with Vi^cj)) := |m^0^ and with initial conditions 
chosen to indicate different types of behaviours. We use geometrized units (k = 1 and speed 
of light equal to 1). Introducing an arbitrary length scale a, various quantities has following 
dimensions. 

Hm ~ a , p</, ~ , ~ a° , m ~ d^-^ , t ~ a . (26) 

Scaling the quantities by the appropriate power of a all equations are rendered dimen- 
sionless and integrated numerically. Since we use the usual classical form of the matter 
Hamiltonian without the corrections from the inverse volume, the length scale d is suitably 
large (eg d > 100£p). 

For the plots shown below in figure(^, figureQ and figure©, the dimensionless scale 
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factor assumed to be 100 and the dimensionless mass is taken to be 0.001. Three dif- 
ferent values of the discreteness parameter are chosen: e = 0.0, 0.5, 0.7. Three differ- 
ent initial conditions for the scalar field and its momentum are taken namely {4>,p^) = 
(1, 10000), (10, 10000), (50,50). As can be expected, the initial value of the scalar field has 
stronger effect on the evolution. For smaller value of the scalar field, the effect of non-zero e 
is virtually absent while for larger values one sees the multiple re-collapse/bounce possibili- 
ties. For the plots, spatially flat model is considered (77 = 0). Only the evolution of the scale 
factor is shown in the figures. The scalar field typically increases first and then decreases to 
smaller values. Both the increase and the decrease of the scalar field are steeper for larger 
values of e. The scalar field momentum also shows similar behaviour. The decrease however 
is much sharper and occurs at later times for larger values of e. 

y I 1 1 1 1 1 1 1 

£ = 0.0 

£ = 0.5 




01 2345678 

log(time) 

FIG. 1: The plot is for the initial values = l,p0 = 10000. The discreteness corrections are very 
small. 
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D) 

^ 6 
5 
4 

01 2345678 

log(time) 

FIG. 2: This plot is for tlic initial values (p = 10,^0 = 10000. For late times, deviations due to 
non-zero e are clearly visible. 

B. Phenomenological Matter 

In fact, it is rather easy to consider more general case of phenomenological matter if one 
proceeds somewhat naively. By this we mean the following. 

In the current treatment, the matter is also supposed to be described by a Hamiltonian. 
However, for the usual 'dust' and 'radiation' which are supposed to represent non-relativistic 
and relativistic matter respectively, one does not write a Hamiltonian. The underlying 
matter dynamics (with time scale shorter than that of expansion of the universe) is supposed 
to ensure a thermal equilibrium. Such matter is then thermodynamically described in terms 
of an energy density, p and pressure, P, with the equilibrium property implying an equation 
of state, P — cup. The first law of thermodynamics in conjunction with adiabatic expansion 
implies the 'continuity equation' — — 3p(l + to) which can be integrated to give, 

p{a) = pe-^^-a^^^^^"''^^^ = p(^^^ ^^^^'^^ , when u; is constant. (27) 
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FIG. 3: The initial values here are: (f) = hQ^p^j) = 50. For late times the evolutions show multiple 
bounces and re-collapses for non-zero e. 

The equilibrium property also implies that the density and pressure are homogeneous. 
Such matter in (approximate) equilibrium is coupled to gravity via a stress tensor of the 
form of the stress tensor of a prefect fluid (and the matter is correspondingly referred to 
as a perfect fluid) which is of course the most general form of stress tensor in the context 
of homogeneity and isotropy. This stress tensor is conserved by virtue of the continuity 
equation implied by adiabatic expansion and the flrst law of thermodynamics. This way of 
coupling matter to gravity is thus consistent with the Einstein equation. 

It is conceivable that some of the mechanisms (eg microscopic matter dynamics) involved 
in establishing thermal equilibrium may be modified, especially if the expansion time scale 
becomes comparable to that of the matter processes responsible for estabhshing thermal 
equilibrium. However, for the large volume regime we are considering, one may justifiably 
assume that the microscopic matter dynamics are unchanged and lead to thermal equilib- 
rium with the usual equations of states. If this assumption is granted then such sources 
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3 

can be incorporated by taking Hm{a) '■= p(a)^, with the scale factor dependence of matter 
Hamiltonian exphcitly specified via that of the t her mo dynamical energy density, p{a). In 
view of the expressions in (|^. this is equivalent to putting Hm = Pcs^- The assumption 
made above amounts to viewing the effective equation of state, := Pes/ PeS being deter- 
mined from the usual equation of state, u. The effective density and pressure are thus not 
ascribed any thermodynamic origin, but are viewed as prescribing a modified coupling of 
the matter to gravity. 

The evolution of the scale factor is now obtained by solving just the (modified) Friedmann 
equation for various choices of a; (a). For example, for constant equation of state, cu, 
and for expanding universe. 



+ 



upa^ 
6 



a\ -i-3w 



Kpa? 



(28) 



Introducing dimensionless quantities. 



a 



upo? 



6 



a 



T 



t 

a 



the Friedmann equation can be written as, 



d 

dr 



(flat model) 
(closed model) 



(29) 



It is convenient to parameterize the equation of state variable as, uj{n) :- 



3— n 
3{n-l) 



(30) 
(31) 

so that 



n = corresponds to the cosmological constant, n = 2 corresponds to the radiation while 
n = 3 corresponds to the dust. 

For the fiat models, the substitution ^ := (^coshA)""^, leads to. 



a /2\"- 



n — 1 \ ere 



(r - To) 



;coshA)"rfA 



2-" J2 "Cr 



r=0 



sinh(2r — n)A 

2r — n 



(32) 



For even integer n, the parenthesis in the r = n/2 term in the summation, is to be 
replaced by A. We have thus obtained the solution of the Friedmann equation (l3Uj) in a 
parametric form, in particular for the cases of interest namely n = 0,2, 3. 
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To recover the e — > solutions, one has only to note that for non-zero n, the left hand 
side of (jH^ diverges for t > tq (say), and implies that one must have A very large. The 
leading term on the right hand side is then ^^e"'^. Likewise, for large A, ^ ~ (^^)"~^- One 
can now eliminate ee^ to get a(t) oc (t — to) " • Clearly, these match with the usual (e = 0) 
solutions. The case of the cosmological constant (n = 0) is discussed in the next subsection. 

For the closed models, again one can obtain parametric form of the solutions which 
involves an integration. However the qualitative behavior can be seen easily. Now a different 
substitution is convenient. Putting ^ := {-^^^)^"'^^^ one obtains, 

r-A 



n — 1 



-a 



A-n 



(r(A) - T„) 



r\y(c^yr^ (33) 

Jxo Jl - (£) (sinhA')2 



Clearly, A must be bounded to ensure < sinh^A < ^. One can also check that the 
integrand is integrable at the maximum value of A. Taking tq = = Aq and r := t^^x 
for the upper limit of integration to be A^ax (defined by sinh^Amax = ^), one obtains an 
oscillatory universe with a finite period given by T := 2rinax for all the three cases of interest. 
The e — * limit is simpler in this case since ^ does not depend explicitly on e. Also Amax 
diverges and therefore for any finite A we can simply take e = in the integrand. For 
n = 2, 3, then one recovers the usual (parametric) form of the solutions. 

For dust and radiation, the corrections due to non-zero e are extremely small. This is to 
be expected since the densities decrease with the scale factor, making the a~^H^ term very 
small. 



C. Cosmological Constant 

This is obtained by taking Hm '■= Ap^^"^ = and as mentioned before, mathemati- 

cally it corresponds to n = 0. In the equation (j^ . we replace p by A. 

For e = 0, one has the well known solutions: a{t) = a(0) exp ( a/ A /6 t) for t] = and 
a{t) = ye/A cosh(v/A76 (t - to) ) for r/ = 1. 

For e 7^ 0, the solutions are qualitatively different. For the spatially flat case, rj = 0, one 
has a re-collapsing solution given explicitly as, 

n 24 

/,\ "max 2 /oA\ 
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The scale factor vanishes only for t ±00 and the solution is non-singular. As e — 0, the 
maximum value for scale factor diverges and one is either in the expanding or contracting 

phase. One can obtain the e = solution by taking the limit e — 0, to 00 with eev 6 
held constant. 

For the closed model, t] = 1, the scale factor is bounded from both above and below as 
1 < < (1 + and the universe keeps oscillating between these values in finite 

time, with the period of oscillation given by, 

T = 2./6/A / , , sinh(A^ax) := 2/e . (35) 

-^0 -^1 - f sinh^A 

As e — > 0, the maximum value of ^ diverges, the period diverges. The usual de Sitter 
solution is recovered by just taking e = in the integrand of equation (jHHj) which gives A ~ 
cr{T — To) and eliminating A. In both the cases above, one can check from the Raychoudhuri 
equation that smaller scale factor is indeed a local minimum while the larger one is a local 
maximum (as a function of t). 

For the cosmological constant as the only source, for e = 0, the equation of state, w : = 
P/p = — 1 remains independent of the scale factor and the Raychoudhuri equation then 
implies an always accelerating universe, precluding the possibility of re-collapse. With the 
e corrections incorporated, the effective density and pressure is modified and the modified 
equation of state acquires a scale factor dependence which is such that there is always a 
decelerating phase before re-collapse. 

Had we taken e — 0, we would have got the classical GR Hamiltonian leading to the usual 
exponentially expanding universe without a re-collapse {1] = 0). The difference equation on 
the other hand shows that discrete wave function, in any one sector begins to oscillate over 
the quantum geometry scale ~ f/^oT-^p for large values of the triad, implying a breakdown 
of pre-classicality (no small scale variations) of the exact solutions. While deviations from 
pre-classicality are expected in the small volume regime due to quantum effects, deviations 



arge volume regime have no such physical reasons and are viewed as an 'infra-red 



in the very 

problem' all^- ^ typical example of such large volume deviations is that of cosmological 
constant (and non-zero spatial curvature). No quantity which has a local interpretation 
such as energy density (A) is ill-behaved and the integrated quantity diverges due to inte- 
gration over infinite volume. In such cases, the homogeneous idealization (integrand being a 
constant) is thought to become in-applicable beyond a finite volume. Noting that the wave 
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function has an oscillation length of ( Aa) ^ js' , one can obtain a bound on the largest scale 
factor, a, by requiring the oscillation length to be larger than the scale of slow variation, 



In the effective picture, we have a maximum value, amax of the scale factor due to re- 
collapse and for ?7 = it is given by a^^^ = in equation (jHH). Demanding that Omax < o 
gives a bound on the cosmological constant: Aip < e/8. For such a cosmological constant, 
the effective picture can be consistent with pre-classical behaviour of exact solutions. Inci- 
dentally, for the currently favored value of the cosmological constant, one has Aip ~ 10~^^°. 

We would also like to note that in [3| re-collapse of a A— dominated universe has been 
obtained from a different perspective. These authors note that the discrete wave functions 
cease to show oscillatory behaviour beyond a critical volume. The value of the critical 
volume matches exactly with the Omax above thus corroborating the WKB intuition ^. 

In the above discussion, we have used the large scale factor expressions for the equations, 
mainly because we were interested in seeing the re-collapse possibility. Secondly we also used 
the effective classical picture right close to the WKB turning point. For inferring qualitative 
implications as done above, this is okay. The full effective classical equations together with 
the domain of validity of the WKB approximation have been given in the previous section. 
The main lesson is that the corrections due to discreteness can change some of the solutions 
qualitatively while for others these are perturbative in nature. 

V. DISCUSSION AND SUMMARY 

Let us recapitulate the key steps in the derivation of the effective Hamiltonian and clarify 
the role of WKB approximation and its interpretation. 

We began from the quantum Hamiltonian constraint While it has a variety of solu- 
tions, we assumed that there exist a class of solutions which are Taylor expandable, at least 
over some bounded interval(s) of the triad eigenvalues. This allowed us to introduce the 
6C±,6^± in equations (jHl 1^ ITUj). Next, we made the further assumption of slowly varying 
amplitude and phase to truncate the 6C±, 6^± to first derivative of the amplitude and second 
derivative of the phase. Now the WKB approximation is made to arrange the expansions of 

^ GD would like to thank Alejandro Perez and Kavin Vandersloot for e-mail exchanges regarding (ist- 
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the trigonometric functions according to powers of h to arrive at a Hamilton- Jacobi equation 
(fT^ for the phase, from the (o(^°)) terms as well as an equation for the amplitude (fTSjl from 
the {o{h)) terms. The self-consistency of the slow variation of the amplitude and the phase 
can be used to infer the conditions for validity of WKB approximation. These were further 
simplified to equations (|15l ITHj) in the region of interest {p ^ q). The effective Hamiltonian 
is read-off from the simplified Hamilton- Jacobi equation ()15|). We also note that the usual, 
second order Wheeler-DeWitt differential equation has been by-passed. 

Since the higher derivatives of the amplitude and phase are also higher orders in h, the 
WKB approximation (keeping leading powers of h) is correlated with the assumption of 
slow variation. As is well known, the WKB approximation fails near the 'turning points' 
at which the first derivative of the phase vanishes. For definiteness, let us assume that 
these occur at some pmin < Pmax- The domain of validity of the effective Hamiltonian is 
necessarily within these values. Since the WKB solutions themselves are only approximate 
solutions of the fundamental difference equation, it is possible that these solutions fail to 
be a good approximation to the exact solutions. This can happen, for example, when all 
the exact solutions fail to be slowly varying (while WKB solutions are slowly varying). Let 
us assume that some particular exact solution is well approximated by the WKB solution 
in an interval (Pmin)Pmax)- The domain of relevance of the effective Hamiltonian is then 
the intersection of these two intervals. The intersection must be non-empty since there do 
exist exact solutions which are slowly varying exhibiting classical behavior. Clearly, the 
largest possible domain for effective Hamiltonian is the interval between the two turning 
points: Pmin,Pmax- Siuce, as explained in the section (jH]), exact solutions are countable 
sums of independent solutions of the difference equation, it should be possible to choose 
linear combinations which approximate any given WKB solution over the WKB domain, 
for example as a least square fit. Precise and useful quantification of a WKB solution 
approximating an exact solution however needs a separate work. 

To construct WKB solutions, one appeals to the connection of Hamilton- Jacobi equa- 
tion with an associated (classical) Hamiltonian dynamics. It is well known (in a gen- 
eral context) that any given solution, $(q'*), of a Hamilton- Jacobi equation selects a sub- 
class of Hamiltonian trajectories in the phase space, namely, those with initial conditions 
(g*(0),pj(0) := |p-|o)- Conversely, a local solution, in a neighbourhood of a point (^q), is 
constructed as $( (?*(t, '?o) ) •= /o l\'t'),q\t') )dt' . Here, q^{t^ql^q}y) denotes a family 
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of solutions of the Euler-Lagrange equation (obtained from the solutions of the Hamilton's 
equations of motion), with initial conditions ql),qo- Thus, the solutions of the effective 
Hamiltonian equations of motion can be used to construct corresponding (local) solutions 
of the Hamilton- Jacobi equation and hence the amplitude and the phase. In this manner, 
the effective Hamiltonian can be used to construct solutions which are to be approximated 
by exact solutions of the fundamental equation ©. 

In summary, at least in the context of isotropic LQC, the discreteness corrections are 
manifested at the level of the effective Hamiltonian. These corrections arise due to a dis- 
crete geometry underlying the continuous one and are controlled by the parameter /iq and 
the Barbero-Immirzi parameter 7. For the cases where a~^Hm grows with the scale fac- 
tor, the corrections alter the qualitative behaviour while in other cases the corrections are 
perturbative in nature. 
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